Nodal/ Antinodal Dichotomy and the Two Gaps 
of a Superconducting Doped Mott Insulator 
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We study the superconducting state of the hole-doped two-dimensional Hubbard model using 
Cellular Dynamical Mean Field Theory, with the Lanczos method as impurity solver. In the under- 
doped regime, we find a natural decomposition of the one-particle (photoemission) energy-gap into 
two components. The gap in the nodal regions, stemming from the anomalous self-energy, decreases 
with decreasing doping. The antinodal gap has an additional contribution from the normal com- 
ponent of the self-energy, inherited from the normal-state pseudogap, and it increases as the Mott 
insulating phase is approached. 

PACS numbers: 71.10.-w,71.10.Fd,74.20.-z,74.72.-h 



Superconductivity in strongly correlated materials 
such as the high-T c cuprates has been the subject of in- 
tensive research for more than twenty years (for a re- 
view see, e.g.,[l|). From the theoretical side, low-energy 
descriptions in terms of quasiparticles interacting with 
bosonic modes have been widely studied starting from 
the weak correlation limit. A different approach views 
the essence of the high-T c phenomenon as deriving from 
doping with holes a Mott insulator 0] ■ The strong corre- 
lation viewpoint has not been yet developed into a fully 
quantitative theory and whether the weak- and strong- 
coupling pictures are qualitatively or only quantitatively 
different is an important open issue. 

The development of Dynamical Mean Field Theory 
(DMFT) and its cluster extensions Q provides a new 
path to investigate strongly correlated systems. These 
methods construct a mean-field theory for Hubbard- 
like models using a cluster of sites embedded in a self- 
consistent bath. Extensive investigations have been car- 
ried out for intermediate interaction-strength using the 
Dynamical Cluster Approximation on large clusters 
The strong coupling limit is more difficult, as only small 
clusters can be employed. Many groups however have 
identified interesting phenomena, such as the competi- 
tion between superconductivity and antiferromagnetism 
the presence of a pseudogap (PG) @, the forma- 
tion of Fermi arcs @, [a j| [lfjf and the existence of an 
avoided critical point In this work we use Cellu- 

lar DMFT (CDMFT) to explore the energy gap in the 
one-particle spectra of the superconducting state when 
correlations are strong. The goal is to identify qualita- 
tive aspects of the approach to the Mott transition in the 
light of recent experimental studies on superconducting 
under-doped cuprates H, 13, 14, TH, lf|, which report 



the presence of two distinct energy scales. 

We consider the two-dimensional Hubbard Model: 



U J2 



(1) 



= etc.;, 



Cj jCr destroys a cr-spin electron on site i, 
is the number operator and ta = fi is the chemical po- 
tential. Only next-neighbor t and nearest-next-neighbor 
t' = — 0.3t hoppings are considered. The on-site repulsion 
is set U = 12t. We implement CDMFT on a 2x2 pla- 
quette. Though this is the minimal configuration allow- 
ing to study a d-wave superconducting state, it already 
presents a rich physics and we think that its deep under- 
standing is an essential step to be accomplished before 
challenging bigger clusters (hardly accessible by the com- 
putational methods presently available). Ti. is mapped 
onto a 2x2-cluster Anderson impurity model which is 
solved using the Lanczos method 0. The CDMFT self- 
consistency condition 18[ is then enforced via the Dyson 
relations £(iw n ) = Q {iu) n ) — G [T](iui n ), which de- 
termines the cluster self-energy S. The hat denotes 8x8 
matrices with cluster-site indices containing both normal 
and anomalous components (Nambu notation). Q is the 
"Weiss field" describing the bath, G[E] is the one-particle 
cluster Green's function [l8[ and tu n = (2n + 1)tt/0 the 
Matsubara frequencies, with j3t = 300. The bath is de- 
scribed by 8 energy levels determined through a fit on 
the Matsubara axis (0 < uj n < 2U), which weights more 
the low frequencies [8j. 

Our main result is the presence of two energy-scales 
on the under-doped side of the phase diagram. We first 
show that this can be established directly from an anal- 
ysis of quantities inside the 2x2 cluster, which are the 
output of the CDMFT procedure. In the left panel of 
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FIG. 1: (Color online). Left: ReEfg vs. uj n . In the inset, 
the U) n — > value as a function of doping 8. Right: The 
distance from the Fermi level (u = 0) of the left (circle) and 
right (square) edge-peaks in the local DOS -ilmfti (see 
inset) are displayed as a function of 8. The dashed line is the 
average of the left and right values. In the inset Gn is shown 
for 8 = 0.06 using a broadening r\ ~ 7 x W~ 3 t to display poles. 



Fig. [T] we display the real part of the anomalous clus- 
ter self-energy £ ano on the Matsubara axis. Only the 
nearest- neighbor component ReE al 2 10 (iw) is appreciably 
non-zero. The main observation is that at low energy 
Ef2°(0) presents a non-monotonic behavior with dop- 
ing 6, as emphasized in the inset. A first characteris- 
tic energy-scale, measuring the superconducting contri- 
bution to the one-particle energy-gap, can be defined 
as Z noc i^fi° '(0) , where Z no( i is the quasiparticle spec- 
tral weight at the nodal fc-points, where quasiparticles 
are well defined. As shown below, and as physically 
expected, Z no d decreases as the doping is reduced to- 
wards the Mott insulator. Hence, Z no dY>^°{Q) decreases 
too due to the behavior of both Z noc i and E al 2 lo (0). We 
stress the sharp contrast of this result with resonating 
valence bond mean-field (RVB-MF) theories [T3|, where 
Z„ O( iI]j2 O (0) corresponds to the spinon pairing amplitude 
which is largest close to half-filling. In the right panel of 
Fig. [T] we show that there is actually another energy- 
scale, which increases when the doping level is reduced. 
This is revealed by looking at the local density of state 
(DOS) -ilmGu in &[£]. In the Lanczos-CDMFT the 
spectral function on the real axis is obtained as a discrete 
set of poles (shown in the inset), which are displayed by 
adding a small imaginary broadening irj. We extract rel- 
evant energy scales by measuring the distance from the 
Fermi level of the gap edge-peaks. While for S > 0.08 
the spectrum is symmetric, an asymmetry appears for 
5 < 0.08. The total energy gap (dashed line in Fig. [J) 
grows with decreasing doping 6, as in RVB-MF theories. 

In order to make contact with experimental observ- 
ables it is useful to obtain momentum-resolved quantities 
from the local cluster quantities. For this we need a pe- 
riodization procedure restoring the translational invari- 
ance of the lattice. Several schemes have been proposed 
Q . Building on previous normal-state studies @, E| we 



use a mixed scheme which is able to reconstruct the local 
cluster Green's function (upon integrating over fc the lat- 
tice Green's function) in the nodal and antinodal points 
better than uniform periodization schemes. Our method 
is based on the idea that, when the self-energies are reg- 
ular the, most suitable choice is to periodize the cluster 
self-energy via the formula 



Y}„(k, lo) = — 



fj, v 



- J /J.I>',(7 



(2) 



(/i, v label cluster sites). The anomalous self-energy 
vjano an( j norma i self-energy E nor in the nodal re- 
gions, where we expect to find quasiparticles, are well 
behaved quantities, therefore we extract them through 
formula l[2"]). In particular, the anomalous self-energy 



acquires 
Eg°(w) (cosfc x 



a standard d^2_ y 2-wave form: E ano (fc,ti;) 



cos k y ). On the other hand, when the 
self-energies develop singularities, the cluster self-energy 
is not a good quantity to be periodized. In Ref. [§], it 
has been shown that this takes place in the normal self- 
energy E nor in the antinodal regions, when the system 
approaches the Mott insulator. In this case, a more suit- 
able quantity to be periodized is the the irreducible two- 

r i ~ a 

point cluster cumulant M^ m {uj) = (w + fi)l — E£ or , 
which is a more local and regular quantity. In the antin- 
odal region, therefore, we can apply formula © to M noT , 
to obtain Ai nor (k,uj) and finally extract the normal lat- 
tice self-energy E nor (fc,u) = w + fx - l/M nm (k,Lu). The 
fc-dependent Green's function can be written as a matrix 
in Nambu's space. 



t k 



CO ■ 



-E ano (fc,w) 



(3) 

The imaginary part of the diagonal entry yields the spec- 
tral function A(k,ui) measured in photoemission. 

In order to compare our results with experiments, it 
is useful to disentangle the normal and superconducting 
contributions to the spectral gap. To this end, we first 
set E ano = in Eq. ©. The results are shown in Fig. 
[2] The fc-points along the nodal and antinodal directions 
are chosen as those where the hig hest peak is observed in 
A(k, w), as done, e.g., in Ref. [lj|. Their actual values are 
shown in the inset of panel C of Fig. [21 Near the nodal 
point (panel A) a quasiparticle peak is well defined at the 
Fermi level (oj = 0) and decreases by decreasing doping. 
In the antinodal region (panel B), a quasiparticle peak is 
also found at the Fermi level for S > 0.08. For 8 < 0.08, 
however, the spectral weight shifts to negative energies 
signaling the opening of a PG, whose size increases as 
(5^0. The behaviour of the PG in the superconducting 
solution smoothly connects to the PG previously found 
in the normal state CDMFT study [§] . The approach to 
the Mott transition is characterized by a strong reduction 
in the area of the nodal spectral peak Z no d, which is plot- 
ted in panel C (green circles). We also plot the area of 
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FIG. 2: (Color online) Spectral function A(k, ui) for different 
8. Broadening rj = 0.03i. Panel A: nodal quasiparticle peak; 
Panel B, normal component (set E ano = in Eq. ©) of the 
antinodal quasiparticle peaks; Panel C, nodal and antinodal 
quasiparticle weights. The inset shows the ^-positions of the 
nodal and antinodal points; Panel D, spectra at the antinodes. 



the antinodal peak Z anoc [, which shows a constant value 
upon the opening of the PG (5 > 0.08). In panel D, we 
restore £ ano , and examine the actual superconducting so- 
lution. The superconducting gap opens in the antinodal 
region (the nodal region is practically unaffected). For 
S > 0.08 the spectra are almost symmetric around the 
Fermi level, as in a standard BCS d-wave superconduc- 
tor. In contrast, close to the Mott transition the PG, 
which originates from the normal component, is super- 
imposed to the superconducting gap, resulting in asym- 
metric spectra. This reveals the origin of the left/right 
asymmetry in the cluster DOS discussed in Fig. [TJ 

In the nodal region the quasiparticle peaks are 
well defined at all dopings and we can expand the 
self-energies at low frequencies. The quasiparti- 
cle residue (1 — 9 w ReEfc(w)) _1 |^_ Q (blue crosses in 
panel C of Fig. [2J numerically coincides with the 
area of the quasiparticle peak Z no d- From Eq. 
|3]), we get A(k,w) ~ Z nod S (w - ^v 2 nod k\ +v 2 A k 
where v nod = Z nod \\J k (t k - E nor (fc, 0)) | and v A 



V^^norfSf 2° (0) | sin k no d\ are the normal and anomalous 
velocities respectively perpendicular and parallel to the 
Fermi surface. v\ physically expresses the superconduct- 
ing energy-scale discussed in the left panel of Fig. [TJ 
We display them as a function of doping 5 on the left 
side of Fig. [3] v no d does not show a special trend for 
5 — > and it stays finite, consistently with experiments 
[20l | . The anomalous velocity, «a *C v nod presents a 
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FIG. 3: (Color online) Left: v no d and «a as a function of 
doping 8 (a is the lattice spacing). Right: low- frequency co- 
efficients of local DOS 7 and of the Raman B2 9 and superfiuid 
density response a, renormalized by the value at 8 — 0.16. 



dome-like shape. This behavior (confirmed by continuous 
time quantum Monte Carlo (CTQMC) calculations (2l| ) 
is in agreement with recent experiments on under-doped 
cuprates showing that, contrary to the antinodal gap, the 
nodal gap decreases by reducing doping [12, [H, 13 • 

The low-energy behaviour of several physical observ- 
ables in the superconducting state is controlled by nodal- 
quasiparticle properties and hence can be related to 
Vnod, va and Z no d- Two specific ratios are particu- 
larly significant, namely: 7 = Z nod / {v no dV a) and a — 
Znodl( v nodV&)- The first one is associated with the low- 
energy behaviour of the local DOS measured in tun- 
neling experiments: N(ui) — A(k, lo) ~ jlu (for 
lu — > 0). Neglecting vertex corrections [12[, the sec- 
ond one determines the low-energy Big Raman response 
function x"( w ) °c auj and the low-temperature (T — > 0) 
behaviour of the penetration depth (superfiuid density) 
p s {T) — p s {0) oc aT. We display a and 7 in the right 
panel of Fig. [3] as a function of 5. a is monotonically 
decreasing (see also CTQMC results [2ljj ) and, on the 
under-doped side 5 < 0.08, it saturates to a constant 
value, in agreement with Raman spectroscopy [l2| and 
penetration depth measurements [22|. Also 7 neatly de- 
creases in going from the over-doped to the under-doped 
side, but it presents a weak upturn for low doping. The 
low- frequency linear behavior of N(li) is well established 
in scanning tunneling experiments on the cuprates [23| . 
However, it is not currently possible to determine the ab- 
solute values of the tunneling slope a from experiments, 
hence the behavior we find is a theoretical prediction. 

We finally turn to the one-electron spectra in the antin- 
odal region, shown in Fig. 31 physically interpreting the 
cluster energy-scales observed in Fig. [TJ We evaluate 
the antinodal gap in the superconducting state A to t by 
measuring the distance from the Fermi level (ui = 0) at 
which spectral peaks are located (panel D of Fig. [2]) . A tot 
monotonically increases by reducing doping, as observed 
in experiments. The data of panel B in Fig. [51 where 
vj ano = 0, allow us to extract the normal contribution 
A„ or . We notice that the peaks found there at negative 
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FIG. 4: (Color online) Antinodal energy gap A to t (circles), 
obtained from the spectra of panel D in Fig. [2] as a func- 
tion of doping S, and decomposed in a normal contribution 
A nor (squares), obtained from panel B in Fig. [2] and in a 
superconducting contribution A sc (diamonds). 



ground-state, or it requires some further ingredient to 
be stabilized against competing instabilities (above all 
antiferromagnetism at low doping ||) remains an impor- 
tant open question addressed to future developments. We 
think however that the assumption of a d-wave super- 
conducting ground-state is a reasonable starting point, 
and the importance of our 2 x 2-plaquette-CDMFT result 
stands in the natural explanation it provides of the prop- 
erties of under-doped cuprates. 
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and M. Le Tacon for useful discussions. M.Ca. was sup- 
ported by MIUR PRIN05 Prot. 200522492, G.K. by the 
NSF under Grant No. DMR 0528969. 



frequency ui pg do not represent Landau quasiparticles in a 
strict sense, but we can estimate the PG as \ui pg \. We also 
display the anomalous contribution to the antinodal gap 
A sc = y/ Af ot — A^ or , and find that, within numerical 
precision, A sc ~ Z anod \E ano (k anod , u) pg )\. The appear- 
ance of A nor signs a downturn in A sc . We interpret A tot 
as the monotonically increasing antinodal gap observed 
in cuprate superconductors, while the superconducting 
gap A sc , detectable as the nodal-slope va (Fig. [3]), is 
decreasing in approaching the Mott transition. 

The concept of two energy gaps with distinct doping 
dependence in the cuprates has recently been brought 
into focus from an analysis of Rama n sp ectroscopy [12j , 
and photoemission experiments [lH , Il4| . which have re- 
vived experimental and theoretical debate [l6| . Our the- 
oretical dynamical mean-field study of superconductivity 
near the Mott transition establishes the remarkable co- 
existence of a superconducting gap, stemming from the 
anomalous self-energy, with a PG stemming from the 
normal self-energy. This is reminiscent of slave-boson 
RVB-MF of the t — J model [H [H, which uses or- 
der parameters defined on a link and includes the pos- 
sibility of pairing in both the particle-particle and the 
particle-hole channels. Compared to the self-energy of 
the RVB-MF, the CDMFT lattice-self-energy has con- 
siderably stronger variations on the Fermi surface [9j and 
additional frequency dependence, which makes the elec- 
tron states near the antinodes very incoherent even in 
the superconducting state. Furthermore, in the RVB-MF 
theory the anomalous self-energy monotonically increases 
by decreasing doping, in contrast to our CDMFT results 
which reveal a second energy scale associated with su- 
perconductivity, distinct from the PG, which decreases 
with decreasing doping. Whether this feature survives 
in larger clusters, representing a property of the real 
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